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SECTION 1 


INTRODUCTION 

PARS is an efficient optimization program for resizing structures 
subject to stress, displacement and flutter constraints. The optimization 
method and the static sensitivity calculations are developed in Ref. 1. 

The aerodynamic and flutter analysis and corresponding sensivity calcu- 
lations capabilities of PARS are based on the WIDOWAC program [2]. PARS 
is based on the SPAR finite element system [3] and like SPAR it is composed 
of a group of individual processors. These processors are used in a 
logical sequence to perform a desired analysis or synthesis task. The 
processors are able to communicate directly and automatically with a 
"Data complex" which contains the information generated by other pro- 
cessors during a computer run. A general characteristic of the technical 
processors is that they are efficient with respect to both core store 
and central processing time requirement. PARS provides an efficient 
code for sizing large or small scale finite element models in the 
presence of strength, thermal and aeroelastic constraints with minimum 
and maximum bounds on structural dimensions. 

Three standard built in looping systems , provided for carrying 
out static strength optimization, flutter optimization, are available. 

The optimization task can also be achieved independently by calling a 
few processors, examining the intermediate results and then entering 
. into the loop or continue calling each individual processors 
sequentially. 
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The optimization procedure is based on a mathematical programming 
method Using the sequence of unconstrained minimization technique (SUMT) 
in which design requirements are incorporated using on extended - 
interior penalty function approach. Approximate analysis techniques, 
design variable linking, constraint deletion and approximate expressions 
for the second derivatives of the constraints are used in conjunction 
with Newton’s method to reduce computational time and increase the order 
of the structural design problem that can be solved. In addition PARS 

may be used with the stress ratio technioue to obtain a fully stressed design 

when the structure is subiect only to stress constraints. Tuo versions of 
stress ratio technique are available to the user. The first one is the well 

known fully stressed design (FSD) and the second is the thermal fully stressed 
design, TFSD [4] . 

PARS can be used to resize structvtral systems that can be idealized 
using truss, beam, shear panels, and triangular or quadrilateral 
membrane or plate elements. All elements except the beam elements may 
be resized. It is assumed that this manual would be used as a companion 
to the SPAR Reference Manual [ 3 J and that the reader is familiar with 
the use of the SPAR program and its terminology for input/output and 
data complex. 

1-1 PARS processors and their use 

There are several ways in which PARS could be used by the structural 
designer. The following are a few suggested sequences for calling SPAR 
and PARS processors to accomplish analysis , sensitivity and resizing tasks. 

(i) Calculation of derivatives of displacement and stresses with respect 
to one or several design variables. 
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a. Use SPAB to calculate displacement and stresses (processors 

TAB, ... GSF) 

b. @ XQT DESVM - define design veuriables 

c. § XQT DKS - calculate derivatives of stiffness matrices (needed 

only if design variables control plate bending elements) 

d. § XQT DGRAD - derivative of displacements with respect to 

thickness or area variables 

e. @ XQT TGRAD - derivatives of displacements with respect to 

temperat\xres design variables 

f. § XQT S GRAD - derivatives of stresses 

(ii) Flutter Analysis 

a. Use SPAR to calculate vibration modes and frequencies 
(processors TAB, ELD ... EIG) 

b. @ XQT AERO - define aerodynamic model 

c. © XQT SUBKRN - calculate aerodynamic matrices 

d. § XQT GAF - calculate generalized aerodynamic forces 

e. © XQT FLT - flutter analysis 

(iii) Sensitivity analysis for stresses, one resizing move and 
reanalysis 

a. Use SPAR to calcvilate displacement and stresses 

b. © XQT DESVAR 

c. © XQT DKS (if resizing plate bending elements) 

d. © XQT DGRAD 

e. © XQT TGRAD (if using temperature design variables) 

f. © XQT SGRAD 

g. © XQT DCON - calculate displacement ratio for displacement 

constraints, if any 
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h. § XQT SCON - calcvilate stress ratios 

i. 6 XQT NEWTON - one resizing move, new design variables. 

§ XQT KSDM - update finite element model for new design 

variables 

k. Use SPAR processors starting with E, EKS, to calcxiLate new 
displacement and stresses. 

(iv) Flutter analysis, one resizing move and reanalysis 

a. Use SPAR to calculate displacements and stresses 

b. @ XQT DESVAR 

c. § XQT AERO 

d. @ XQT SUBKRN 

e. § XQT GAF 

f. § XQT DKM - calcvilate derivatives of generalized stiffness and 

mass matrices 

g. @ XQT FLT - flutter analysis and derivatives 

h. § XQT NEWTON - one resizing move, new design variables 

i. @ XQT FLT - new flutter analysis, no mode update 
or 

i. @ XQT KSDM - update finite element model for new design variables 

j. Use SPAR to generate new modes (processors E, EKS through EIG) 

k. @ XQT GAF 

l. @ XQT FLT - new flutter analysis, after mode update. 

(iv) Automated static optimization 

a. @ XQT TAB ■ 

b. @ XQT ELD 

c. § XQT TOPO 

d. 6 XQT DESVAR 
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e. § XQT SOL (with IDP = 28, See SOL), (for mathematical programming) 

or , 

e. g XQT FSDL (with IDP = 28, see FSDL), (for PSD or TFSD resizing) 

(v) Automated Flutter Optimization 

a. § XQT TAB 

b. @ XQT ELD 

c. @ XQT TOPO 

d. § XQT DESVAR 

e. @ XQT AERO 

f . @ XQT SUBKRN 

g. @ XQT FOL (with IDP = 28, see FOL) 

1*2 Input Rules 

The input rules used for SPAR apply to the PARS processor with one 
addition that of a namelist type input. The namelist input is very 
similar to RESET input except that it provides also for array input. 

The following are namelist input rules. It is assumed that the namelist 
name is XX. 

The first character in card to be read is always ignored. The second 
column of the first card must be a $ immediately followed by the namelist 
name and a blank. If the name is not XX, the input medium is searched 
for $XX. The remainder of the first card and following cards may 
contain any combination of the legal data items which are separated 
by commas (a comma after the last data group is ignored). The last 
input record is terminated by a blank followed by $END. 

The forms the data items may take are: 



Varisible Name = Constant 


Subscribed Variable = Constant 
Subscripted Variable = Set of Constants 

The set is represented by constants separated by commas, or k* 
constant where k* constant may represent k constants (k is an unsigned 
integer constant). 

Array Name = Set of Constants 

The set is represented by constants separated by commas, or k* 
constant may represent k constants (k is an unsigned integer). 

The data items that appear on the input records need not appear 
in the same order as the corresponding variable or array names in the 
namelist. All variable or array names of the namelist need not have a 
corresponding data item in the input records; if none appears, the 
default values will be Tised. 

Blanks must not be embedded in a constant or a repeat constant \ 
field, but may be used freely elsewhere in a data record. The name of 
an array and the value of its first elements must appear on the same 
record. The last item on each record that contains data items must be 
a constant followed by a comma. The comma is optional in the record 
that contains or precedes the $END sentinal. 

Example : 

Col. 2 

First Data Card: $NAM1 1(2,3) = 5»7»9» 

Second Data Card: 0,10,L(3) = ^t.3, A=lt,3, 

Third Data Card: 8 * U. 3> J=^».2, 
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Last Data Card: 


$END 



If the data is vised with the namelist example the following action takes 
place: the first data card is read and examined to verify that its name 

is the same as the namelist name. The integer, constant 5 is placed in 
1(2, 3), 7 in 1(3,3) and 9 in l(U,3). When the second data card is read, 
integer 0 is placed in 1(5,3), and 10 in 1(6,3); the real constant h.3 
is converted to integer U and placed in L(.3)» integers U and 3 are 
converted to real constants and placed in A(l) and A(2), respectively. 
When the third data card is read, the real constant U.3 is stored in 
A(3), A(U),...,A(10); and the real constant U.2 is converted to integer 
U and placed in J . 
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SECTION 2 


DESIGN VARIABLE DEFINITION AND INTERFACE WITH FINITE ELEMENT MODEL 

Processors DESVAR and KSDM described in this section, establish the 
relationship between design variables and the finite element model and 
provide the means of changing that model as the design variables are 
resized. Defining the design variables is the first task for a PARS user. 
Sometimes there is a simple one-to-one relationship between design 
variables and parameters of the finite element model. For exan 5 >le, the 
design variable may be the thickness of a panel and all the plate finite 
elements are set equal to the value of that design variable. The 
relation could be more complicated, for example if design variable 
linking is used. The PARS program allows for any general linear 
relationship between design variables and the resized parameters of the 
finite element model. These resized parameters may be thickness or 
cross sectional areas of finite elements, nodal temperatures and nodal 
masses. Some elements in SPAR (such as the beam elements and the fluid 
elements) cannot be controlled by design variables. Thus, if such 
elements are used in a model, their structural dimensions remain 
unchanged during optimization. 

The thickness and cross sectional areas of finite elements are 
controlled in SPAR by section property input (SA, SB, etc.). Consequently 
in PARS the design variables control section properties and all the 
elements that have the same section type and nmber have the same 
section properties. 
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The following is a list of the types of elements belonging to SPAR 
element library which can be resized, and their section type. 

Name Description 

E23 Bar elements having only axial stifftiess , section type BC. 

E 31 Triangular membrane, constant-stress element, section type SA. 

E32 Triangrilar bending element, section type, SA. 

E33 Triangvilar membrane + bending element, section type SA. 

Ehl Quadrilateral membrane, element section type SA. 

EU2 Quadrilateral bending element, section type SA. 

EU3 Quadrilateral membrane + bending element, section type SA. 

EU1| Quadrilateral shear panel, section type SB. 

The resized dimension for section type BC is the bar cross sectional 
area and for section types SA or SB it is the plate thickness. 

2.1 DESVAR Design Variable Definition 

Function: In a structural optimization problem, it is necessary to define 

the design variables which control various parameters of the structure. 

DESVAR permits the user to relate a structural parameter q as a linear 
function of the design variables in the following form, 

n 

q = n + Z V. F. (2.1— l) 

’ ^ 1=1 ' ^ 

where, 

q is a structural parameter. 

q^ is the minimum value of the parameter 

V. is the value of ith design variable (always positive) 
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is the amount by which a structural parameter can be changed 
by making a unit change in the value of the ith design 
variable, i.e., v^. 

n is the number of design variables* 

Most of the are zerQ, the user has to specify only the nonzero 

ones. 

Also (2.1-2) 

3v^ 1 

A structural parameter, which is controlled hy a design variable can 
he either 

(i) a tuning mass lumped at a grid point or 

(ii) a thickness or an area of a finite element or 

(iii) a nodal temperature at a grid point. 

Parts of the structure may he assigned no design variables and remain 
unchanged during the optimization process. 

Input /Output Data Sets; Input data sets are element and section 
definition sets generated hy TAB and ELD. 

Output data sets are governed hy the input information relating to 
the design variable definitions. 

The following output data sets are written into library DLIB (see 
reset controls); 

1, MASS DV 0 0, contains data relating design variables and tuning 
masses , in the following order nj nv f which indicates that 
design variable nv controls the mass concentrated at joint nj 

by a scaling factor f. NI = 3, NJ = number of controlled masses. 

*At most five nonzero factors are allowed. When PSD or TFSD resizing 
is used each section shoiHd be controlled by a single design vai’iable. 
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2. NTDV SEQ tset 0; contains the identity of the design variables 

that control nodal temperatures. NI ■ 1, NJ “ maximum number of 
variables that control one nodal temperature, 

3. TEMP DERI nset 0, contains the scaling factors which relate 
design variables and nodal temperatures. NI “ maximum number of 
variables that control one nodal temperature, NJ “ number of joints. 

It. NSDV sname 0 ID. where sname is a section name (BC, SA or SB) 
and ID is the section's identification number (ll, 13 or l4) 
contains the number of design variables (max. 5) that control 
each section. NI = number of sections, NJ = 1 

5. KSDV sname 0 ID; contains the identity of the design variable 

that control each section. NI - numoer of sections, NJ = minimum 
number of variables that control a section. 

6. FSDV sname o ID; contains the scaling factors for each section. NI 
number of sections, NJ = maximum number of variables that control 
a section. 

7. DCON BLOCK 0 0 Contains data about element types and design 
variables . 

8. DESl VARl 0 0 Contains the values of the design variables. NI=1, 

NJ = nximber of design variables. 

9. GAGE VARI 0 0 contains the reference values of the design 
variable (often corresponds to the minimtim structural parameter 
controlled by the design variable). NI = 1, NJ = number of 
design variables. 

10- DVMX VARI 0 0 contains the maximum values of the design variables 
NI = 1, NJ = number of design variables. 
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Reset Controls ; 

Name Default Value 

DLIB 1 

SET 999 

BLDV 5000 

BLIB 1 


Description 

Destination Library for design variable 
definition data sets 

The set number associated with design 
variables controlling nodal temperat\rres 
Block length reserved for generating 
design variable definition data sets. 
Library containing data sets generated 
by TAB 


Execution Sequence: 

The following is the sequence of data cards 


eXQT DESV . 
RESET Cards 


Design variable definition cards 
Design Variable Definition: 

The following is the sequence of the design variable definition data cards 
(a) MASS: nj^, nv^, f^: njg. nv^, f^; nj^. nv^, f 3 : • = $ 

The first record identifies tuning MASS parameters. The subsequent 
records form triads that are separated by a terminator mark (: on UNIVAC; 
on CDC). A triad n j , nv, f indicates that the design variable nv 
controls the mass concentrated at node nj by a scaling factor f. 
Example: MASS: 51 » 2, 3«25: 1^» 1.15$ 

25, 3, 2.15: 17, 2, 1.75 $• 

The above data indicates that design variable 2 controls two masses m^^^ 
and m^.^ lumped at nodes 51 17 respectively; design variable 3 

controls the mass design variable U controls the mass ra^j^. Also, 
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(b) SECTNAME: tiv^ isec^, jsec^^j 

nVg, isec^, Jsec^, f^: nv^, isec^, Jsec^, 

nvj^, isecj^, jsec^, fj^: 

The first record in each set of cards contains the section type. The 
subsequent records form quartets each separated by a terminator mark 
(:on UNIVAC; on CDC). A quartet nv, isec, jsec, f indicates that the 

thickness or area of elements controlled by sections isec to jsec 

(assigned in loop limit format, jsec can be omitted if it is equal to isec) 
is controlled by design variable nv with a scaling factor f. Hote 
that sections that are not resized should have higher numbers than sections 
which are resized. 

Example 1 

SA: 3, 2, 5, 1.5: 2, 1, 1.7 $ 

The above data indicates that the thickness t^ of sections 2 to 5 
of type SA is controlled by design variable number 3. Also the 
thickness, tg of section 1 of the same section type is controlled 
by design variable number 2. Also 
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(c) TEMPERATURES: nj^, mj^, nv^, f^ $ 

njg, mjg, nvg, nj^* nv^, f3: $ 

The first record identifies nodal t emper at \ire parameters . A quartet 
nj, mj, nv, f indicates that the nodal temperatures specified at the 
^ioints nj to mj (assigned in loop limit format) are controlled by design 
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variable nv vith a scaling factor f. 


Exan^le 

TEMPERATURES $ 

5, T, 7, 1.175: h6, U6, 8, 5-25: 50, 50, 7, 3.25 $ 

33, 33, 7, 1.00 $ 

The above data indicates that design variable 7 controls the 
temperatures T^, T^, T^, T^^ and T^^ Joints 5 to 7, 50 and 
33 respectively; design variable 8 controls the temperature Tj^g 
with a factor 5.25. 

Also, 


*’■33 

3^=5.25, 3.25, 1.0 and 

o 7 I 


= 1.175, for J = 5, 7. 


(d) DESIGN: nv^, mv^, v^: nVg, mvg, Vg $ 


: nv, , mv. 


k- ''k * 


This specifies the values of the design variables. A triad nv, mv, v 
indicates that design variables nv to mv (assigned in a loop limit format 
have the values v. 


Example 


DESIGN $ 

1, U, 1.76: 5, 7, 2.58: 8, 50, 1.69 $ 

The above indicates that design variables 1 to U are equal to 1.76; 
design variables 5 to 7 are equal to 2.58 and design variables 
8 to 50 are each equal to 1.69. 
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(e) GAGE: nv^, mv^, gj_ : nv^, g^: 

The GAGE card specifies the reference values for the design 

vAri flhi es (see discussion in section 5»l)* reconnnended that the 

j-0f0i*ence values he set to the minimuiii gage of the structural parameters 
controlled hy the design variables. A triad nv, mv, g indicates that 
design variables nv to mv (assigned in a loop limit format) have the 
reference value equal to g. 

Example 

GAGE $ 

1, 15, 0.05: 16, 50, .01 $ 

The reference values for the design variables 1 to 15 are .05 and 
for design variables l6.to 50 are equal to .01. 

If one or more of entries (a) to (c) is not present in the 
structural system then the corresponding input cards can be omitted. 

For example! when themnal constraxnts are not consxdered xn design 
procedures, there is no need to specify TEMPERATURE design variables. 

(f) DMASS: nv^, mv^, dm^: nv^, mVg, dmg-. 

The DMASS card specifies the base values for the mass derivatives, 
which may account for mass contributions such as secondary structure 
which are not included in the finite element model- but are proportional 
to the design variables. A triad nv, mv, dm indicates that design 
variables nv to mv have the base mass dm. 

Example 

DMASS $ 

1, 15, 0.001: l6, l6, .005 $ 

The base mass derivative values to design variables 1 to 15 are .001 
and for design variable l6 is .005. 
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(g) DVMX: mv^» mVg, ym^:..^. 

The DVMX card specifies upper bounds on design variables. A triad nv, 
mv, vm indicates that design variables nv to mv have the maxlBum value vm. 

This data is needed only if upper limits on design variables are desired. Note 
that data items (d) , (e), (f ) , and (g) can be skipped and AUS used instead to 
create the corresponding data sets. 

2.2 KSDM ; Structural parameter update Function: KSDM performs one 

or more of the following functions depending upon the values of the 

reset parameters (see RESET controls) 

1. Updates the values of structural dimensions (areas for one 
dimensional elements, thicknesses for 2 -dimensional elements) nodal 
temperatures and tuning masses which are controlled by design variables ^ 

(see processor DESVAR). The dependence of a structural parameter q 
on the design variables is 


where a - is the minimirni value of the parameter q; 

0 

V : is the ith design variable (> 0.) 

1 

ni is number of design variables 

F : is the amount by which a structural parameter q is changed 
i 

due to a unit change in the Value of ith design variable, v.. 

The user has two options for defining the minimum values q^ of structural 
parameters which are cross sectional properties (such as thickness or areas) 
The first option is to use the structure defined by SPAR processor TAB 
as the minimum gage structure, dliis option is selected by setting the 
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reset parameter MING to -1 the first time KSDM is executed. The other 
option is to use the structure defined by TAB as the initial design 
and to have KSDM calculate the minimum gage values from Eq. (2.2-1) 
using the initial values of the design variables. This option is 
selected by setting MING =0. 

2. To generate the minimxmi-gage values of the thicknesses 

t and areas A for individual section of a given section type, 
o o 

say SA. This calculation is made when the structure defined by SPAR 
is the initial design (reset control MING = O). 

3. To evaluate the derivatives of mass with respect to design 
variables. The derivatives of the mass is the sum of base values 
supplied by the user in processor DESVAR and the derivatives of the 
structural mass. The latter (derivatives of structural mass) are 
generated internally in KSDM. The base values of the mass derivatives 
may be used to account for mass contributions such as secondary 
structiare which are not included in the finite element model but are 
proportional to the design variables. For the design variables which 
control nodal temperatures , the base mass derivatives may represent 
the mass penalty (e.g., insulation) associated with temperatiure 
change. It can also be used to represent the cost penalty associated 
with installation and maintenance of any cooling system device to 
reduce temperature. In the latter case eqxii valent mass penalty is 
used for the base mass derivatives. 

Input Data Sets: 

1. The data set DESI VARI 0 0 (NI = 1, NJ = HDV) where NDV 

•is the number of design variables contains the value of the 
design variables on library DLIB. 
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2. 


The data set DMAS DV 0 0(NJ=1, NI=NDV) or library DLIB contains 
the base mass derivative values (if any) which may account for 
mass contributions such as secondary structure which are not 
included in the finite element model but are proportional to 
the design variables. 

The above two data sets have to be generated by DESVAE or AUS the first 
time KSDM is executed. DESI VARI data set is updated by NEWTON or SCON. 

3» Data sets generated by TAB on library BLIB. 

1*. The element packets from library ELIB. 

5. Data sets generated by DESVAE on library DLIB. 

6. Nodal temperature (if any) data set generated by AUS on 
library TLIB. 

Output Data Sets /generated on library DLIB 

1. CDIM namesect 0 ID (NJ=1, NI=NSECT) contains current structural 

parameters values . 

2. GDIM namesect 0 ID (NJ=1, NI=NSECT) contains minimum gage 
structural parameters values . 

3. ISCT element name 0 ID (NJ, NI=NELEMT) contains the 
corresponding section number of each element. 

It. INGS elemt name 0 0 (NJ=5, NI=NDV) contains the elements which 
are controlled by each design variable and by what factor. 

5. MASS DER 0 0 (NI = NDV* NI - 1) contains the derivatives of the mass 
of the structure with respect to the design variables. 

6. MING BW 0 0 (NI =■ 1, NJ = 1) contains the minimum gage structural 

weight. 
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Reset Controls: 


Name 

BLIB 

ELIB 

GAZEBO 

CIZERO 

TLIB 

BLOC 

DLIB 

MING 


Defavilt Value 
1 


1 


l.E-20 


l.E-20 


1 

10 


1 


1 


Description 

Library containing TAB generated 
data sets. 

Library containing the element 
packets . 

Zero-test parameter, (beam area) 

X (shear modul\os). 

Zero-test parameter, beam non-uniform 
torsion constant . 

Source library for nodal temperatures . 
Maximum mmiber of blocks in INGS data 
set. Shoiild be set to zero to save 
core if FSD or TFSD resizing is used. 
Source library for design data sets 
created by DESVAR and destination 
library for data sets created here. 

MING = -1: The original structure 

generated by TAB) is taken as the 
minimum-gage structure . 

MING = 0: The minimum gage thicknesses 

are determined using Eq. (2.2-1) 

MING >0: It is asstimed that minimm 

gage values for structural parameters 
are already generated in the previous 
executions). 
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Reset Controls (Continued) 


Name 

IND 


ITDV 


MASS 


G 


Default Value 
0 



0 


1 . 


Description 

IND=0; when KSDM is executed alone, 

IND=1; when KSDM is in an optimization 
loop. 

ITDV=0 indicates that temperature design 
variable are used. This reset control is 
not required when KSDM is executed as part 
of SOL or FSDL 

MASS=0 indicates that no tuning mass 
design variables are used. 

MASS?^0 indicates that timing mass design 
variables are used. 

Gravitational constant. Used when mass 
properties are desired while the TAB data 
included weight densities. 
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SECTION 3 


STATIC CONSTRAINTS AND THEIR DERIVATIVES 

The following sections describe PARS processors that calcialate 
derivatives of displacements, stresses, and displacement and stress 
constraints. Processors DGRAD and TGRAD calculate the derivatives of the 
displacements with respect to design variables that control structural 
dimensions and nodal temperatvires , respectively. Processor SGRAD 
calcxilates derivatives of the element stresses. Processors DCON and 
SCON are used to define displacement and stress constraints and calculate 
the derivatives of these constraints with respect to design variables. 

3.1 DKS ; Derivative of K Matrices 

Function: DKS Computes the derivatives of K matrix with respect to 

design variables for plate bending elements. This is needed for 
computing the analytical derivatives of the displacement vectors. 

Input data sets: 

Input data sets are EFIL and design variable definition data sets 
generated by processors E, EKS and DESVAR respectively. 

Output Data Set is 

DERK EXY Itype nnodes: This data set contains the derivatives of 

the stiffness matrix for each element in the intrinsic form. It is 
stored in the following form 
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DERK (NDOF, NEL, NDV) 

where NDOF is the number of intrinsic degrees of freedom 
NEL is the number of elements 
NDV is the nxunber of design variables 


Reset 

Controls : 


Name 

Default Value 

Description 

BLIB 

1 

Library containing TAB generated 
data sets 

ELIB 

1 

Library containing the element 
packets 

KSLI 

1 

Library containing the output data 
sets from this processor 

DVLI 

1 

Source library for design variable 
data sets 


3.2 DGRAD - Gradient of Displacement Generator . 

Function: DGRAD is used to c^culate the derivative of displacement 

vectors with respect to design variables . The design variable is a 
structural dimension which controls the thickness or area of one or 
several finite elements of one or more types . 

In calculating the derivatives of the displacements the dependence 
of thermal or dislocation loads on the design variables is accounted for while 
the dependence of inertia loads on the design variable is not accounted for. 

Input /Output Data Sets: 

Input data sets are usually the joint displacement data sets 
resident in QLIB as STAT DISP iset neon or S0URCE=N1, N2, n^, n^. 

The SOURCE statement names a data set resident in QLIB which 
replaces STAT DISP iset neon as the source of joint motion data. 
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The Output data set is: - GRAD DISP i set neon; where iset, neon have 

eonventional SPAR meanings. It is stored in library GDLIB as an array 
of the following dimensions: 

GRAD (JDF * JT, NDV, NCASE) 

where 

JDF: is the number of aetive degrees of freedom. 

JT: is the total number of joints. 

NDV: is the number of design variables. 

NCASE: is the number of load oases in the speoified iset. (The third 

entry for the data set name). The GRAD DISP iset neon data 

set is generated for ioase=ioasel, ioasel+1, ,iease2. 

Reset Controls : 

Desoription 

So-urce library for the element packets , 
Soxirce library for STAT DISP iset neon 
or the SOURCE data set supplied 
alternatively. 

Load Set (iset). 

Library for design variable data 
sets (DESVAR) 

Library for initial strain data sets, 
icasel (if zero, the first load case 
number used in STAT DISP data set , 
if present, is assumed). 
icaseP (if zero, the last load case 
nvimber used in STAT DISP data set, 
if present, is assigned). 


Name 

ELIB 

QLIB 


Default Value 
1 
1 


SET 

DVLI 

SI LI 
LI 


1 

1 

1 

0 


L2 
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Reset Controls (Continued): 


Default Value 
CON 0 

lEA 0 


GDLIB 1 


KILIB 1 

K K 


Description 

Uth word of STAT DISP iset neon. 

Read error control parameter: 

IEA=0: Both LIB and Fin errors are 

displayed, and run aborted. 
IEA=1: LIB error is displayed only. 

IEA=2: Neither LIB nor FIN error 

are displayed. 

Destination library for the gradient 
of displacement data set GRAD DISP 
iset neon. 

Library containing factored stiffness 
matrix. 

Name of the stiffness matrix. 


Execution sequence: 

The following is the sequence of data cards: 


@XQT DGRA 


RESET Cards 

S0URCE=N1, N2, nl, n2 (trailing words mask - filled if omitted). 
Note that nl will override the SET reset control. 


3.3 TGRAD Derivatives of Displacement Vectors with respect to Design 
Variables Controlling Nodal Temperature 
Function: The processor TGRAD computes the derivatives of displacement 

vector with respect to design variables which control the nodal temperatures. 
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Temperatures specified at several joints can /be controlled by the 
same design variable or conversely a single joint temperature can be 
controlled by more than one design variable. Basic definitions and . 
input information for the temperature design variables are provided in 
the section describing processor DESVAR. 

A nodal temperature Tj at the joint j is related to the design 
variables by 


i=n 


T ,= + Z V. »F.. ; j = 1, m 

0 Oj 1 ij 


where 


T : 
°j 


V. : 
1 


F. 

ij 


nodal temperature at joint j 
reference temperature at joint j. 

This usually corresponds to the minimxmi temperature at the node j 
Value of the ith design variable (always positive) 
the amount by which the temperature can be changed by making 

a mit change in the value of the ith design variable. That is 


3T 


= F 
3v. ij 


For those nodes which are not controlled by temperature design variables 
"ij = ° 


n: the number of design variables, 
m: the number of nodes. 


Output Data Sets : 

Output data sets are: 
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1. EQUF FORC tset iblock: ' This data set contains the derivatives 

of equivalent nodal force generated due to the presence of thermal loadings, 
vith respect to design variables controlling 

2. IS EXI tset iblock: This data set contains the derivatives of 

strain produced due to the presence of thermal loadings with respect to 
design variables controlling T.. 

3. GRAD DISP iset icon: This data set is originally created by 

processor DGRAD. TGRAD fills in empty block spaces of this data set with 
the derivatives of displacement vectors with respect to design variables 
controlling nodal temperatvires . They are stored as 

GRAD (JDF, JT, NDV, NCASE) 

where 

JDF: is the number of unconstrained degrees of freedom. 

JT: is the number of joints. 

NDV: is the number of design variables. 

NCASE: is the total number of load cases in a given load set (iset). 

Input Data Sets: 

The following two input data sets generated in processor DESVAR 
are used by TGRAD 

1. TEMP DERI tset 0. Stores the derivatives of nodal temperature 

vector [T] = [T.] of dimension JT with respect to design variables, JT 
J 

being the mmiber of joints. 

2. NTDV SEQ tset 0: Contains the design variable numbers associated 

with the design parameter controlling the nodal temperatures. 
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Reset Controls: 


Name 

Default Value 

Description 

SET 

1 

Load set nxmiber. 

TSET 

999 

Set number associated with derivatives 
of equivalent nodal force data and 
derivative of initial strain data. 

Do not use same number for RESET 
parameters SET and TSET. 

LI 

0 

First load case number. In default it 
takes the lowest case number previously 
generated. 

L2 

0 

Last load case number. In default it 
takes the highest load case number 
previously generated. 

ELIB 

1 

Library containing element packets. 

INLI 

1 

Source library for all applied load 
data (e.g., nodal or element 
temperatures etc. ) . 

FEFLIB 

1 

Destination library for derivatives 
of equivalent nodal force data. 

ISLIB 

1 

Destination library for derivatives 
of initial stiain data. 

ISBL 

896 

Block length for generating the 
derivatives of initial strain data. 

CON 

1 

Constraint case (see discussion for 
processor INV) 
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Reset 

Cont rols : (Continue d ) 


Name 

Default Value 

Description 

GDLIB 

1 

Destination library for derivatives 
of displacement vectors data. 

K 

K 

Name of stiffness matrix. 

KILIB 

1 

Library containing factored stiffness 
matrix. 

DVLI 

1 

Library containing data sets 
created by DESVAR. 


3.1* SGRAD Stress Gradient Generator . 

Function: SGRAD evaluates the derivatives of element stresses at the 

centroid of each element with respect to the design variables. 

Input Data Sets: 

1. GRAD DISP iset icon: Contains the derivatives of displacement 

vector with respect to design variables. 

2. IS EXY tset mask: Contains the derivatives of initial strains 

with respect to design variables controlling nodal temperatiires 

3. NTDV SEQ tset mask: Stores design variable numbers controlling 

nodal temperature. 

Data sets 2 and 3 are reqxrired if temperature design variables are 
present. 

Output Data Sets: 

1. GOST EXT iset icase: This data set contains the derivatives of 

the stress tensor (at the centroid of each element) with respect to 
design variables in the following form: • 

GOST (LSFE, NEL, NDV) 
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where, 

LSFE; is the nvmiber of stress components generated per element. 
HEL: is the number of the elements in a given element type EXY. 
NDV: is the number of the design veiriables. 

Reset Controls: 


Name 

ELIB 

QLIB 

LI 


L2 


SET 

CON 

LREC 

INST 

GDLI 


Default Value 
1 
1 


1 

MASK 

5600 

25 


Description 

Library containing the element packets 
Library containing the derivatives 
of initial strain data sets. 

First load case nxmber. In default 
it assumes the lowest load case number 
for which derivatives of displacement 
vectors are generated. 

Last load case number. In defaiilt 
it assumes the highest load case 
niimber for which the derivatives of 
displacement vector are generated. 

Load set nvmiber. 

Constraint case nvmiber (see INV 
discussion) . 

Block length required for generating 
the derivatives of stress data set. 

The source library for INST EXY ’ 
mask mask data set. 

Source library for derivatives of 
displacement data. 
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Reset 

Cont rols : ( Cont inue d ) 


Name 

Default Value 

Description 

GSLI 

1 

Destination library for derivatives 
of stress data. 

DVLI 

1 

Library containing data sets generated 
by processor. 

TSET 

999 

Load set associated with temperatvire 
design variables (see DESVAR for 
discussion). 


Execution Sequence 
@XQT SGRA 

Reset cards (if any) 

SOURCE = ^2 words MASK - filled if omitted) 

The SOURCE statement names a data set, resident in QLIB which 
replaces GRAD DISP iset icon, as the source of displacement gradient 
data. 

Subsequent control cards designate specific groups, of elements for 
which gradient of stress data with respect to design variables is to be 
computed, as illustrated below: 

E21: 1+,T: 10$ E21 groups U through 7, and group 10. 

E33$ All E33 elements 

EU3: 1:7:10, 12$ EU3, groups. 1,7,10,11,12. 

El*l»: 1,5:7,17$ etc. 

If the control cards above are not given, derivatives of stresses 
will be computed for all elements which are controlled by design 
vairiables. 
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3.5 DCON : Displacement Ratio Generator 

Fvinction; DCON determines the displacement ratios and the derivatives 
of these ratios at the nodal points of the structure where maximum 
displacements are prescribed. The displacement ratio is the absolute 
value of the ratio of a computed displacement component at a joint to 
its allowable value. 

If Q. is the ith displacement ratio, the displacement constraints 
are formulated in the following manner: 

= 1 _ ^ 0; i = 1, NDC (3.5.1) 

NDC is the number of displacement constraints . 

Input data sets are: 

1. STAT DISP iset icon: Contains the displacement vectors stored in 
SYSVEC format. It consists of n blocks, where n is the 
number of cases in the designated load set. Block 1 is the 
solution corresponding to case 1, block 2 corresponds to case 2, 
etc. 

2. GRAD DISP iset icon: Contains the derivatives of displacement 

vectors with respect to design variables . All the NDV blocks 

of case 1 •are stored first, NDV blocks of case 2 are stored next 

and so on, where NDV is the number of design variables. 

Output data sets are: 

3. DISP CON 0 0: Contains the allowable values of displacements. 

U. EADS CON iset icase: Contains the absolute value of the ratio 

of the computed displacement component at a joint to its 
maximum allowable value. This data set is a vector of dimension 
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NDC where NDG is the maximum number of displacement constraints. 
5. GGl CON set coid - contains the constraint functions given by 
Eq. 3.5.1 for those constraints for which £ GLIM (see reset 
controls) (,i.e. most critical constraints) 


6 . DGl con set coid - contains the derivatives of the constraints 
included in the data set GGl CON. Each block in this data 
set contains -fehe derivatives of a single displacement constraint 
with respect to all design variables. 

Reset Controls: 

Name Default Value Description 

Library containing the element 
information packets, 

Library containing the displacement 
data sets. Also destination 
library for EADS COl? data 
set 


EUB 


QLIB 


GDLI 3 


1 


Library containing the derivatives 


SET 

CON 


LI 


1 

1 

1 


of displacement data, Also destination 
library for GGl CON and DGI CON data 
sets , 

Load set number, 

Constraint case number (see 
discussion for INV processor) 

First load ca.se number. 
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Reset Controls (continued); 
N»mp Default Value 
L2 0 


DISC 25 


GLIM 1 


COID 


1 


Description 

Last load case number. In defaxilt it 
assumes the maximum number of load 
cases for which STAT DISP data set 
exists. 

Destination or source library for 
DISP CON mask mask data set. 

Critical constraint cutoff point. 

Only constraints that satisfy ^GLIM 
are included in the data set GGl CON. 
Constraint identification number. 

Used as the last name of data sets 
GGl CON and DGl CON. Should be 
different than any constraint identi- 
fication number for stress, flutter or 
any other constraints generated by 
other processors . 


Execution Sequence: 

The following is the sequences of data cards. 


@XQT DCON 

Reset cards (if any) 

DISCON: n^, adv^: jg, adVg: n^, ^ 2 * : nj^, adVj^$ 

The DISCON card specifies the allowable values for the nodal displacements. 
The entries in between the two record delimiters (:on UNIVAC) indicate 
that nodal displacement at joint n in direction j has an absolute value 
less than or equal to adv. 
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Example 

DISCON: 35, 3, 0.5T; 35, 1, 0.05$ 

25, 3, O.U5$ 

The allowable values for w^^ and at joint 35 are' 0.57 and 0.05 

units respectively. The allowable value for w^^ is 0.U5 unit. 

3.6 SCON: Stress ratio generator 

FUNCTION: SCON determines the equivalent stress ratios and either 

generates their derivatives with respect to the design variables xdien 
mathematical programming optimization is used, or generates resize 
factors when FSD (Fully Stressed Design) or TFSD (Thermal Fully Stressed 
Design [4] is applied. A design variable can control the magnitude of 
nodal temperatures or the structural dimensions of the finite elements. 

The equivalent stress ratio for the different element types are 
calculated as follows t 

(a) Axial bar element E23: 

SR = •^ (3.6-1) 

is the calculated stress in the element which can be either 
compressive or tensile, and 

is its appropriate allowable tensile XT or compressive XC 
depending on the sign of o^. 

(b) For two dimensional elements E31 E32 E33 

EUI EU2 EU3 

the Tsai Wu theory [5] is used, with the square root of the failure 
criterion used as a stress ratio, 


a 

1 
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( 3 . 6 - 2 ) 


SR = + (yT YC^ ^2 ~ 


(Oif (Og)^ 

XT » XC “ YT * YC 




where are the computed values for the normal and shear 

stresses evaluated at the centroid of the element. 

XT and XC are the tensile and compressive allowables in the.X 
direction and YT and YC are the allowables in the Y direction. is 

the allowable shear stress. The compressive allowables (XC & YC) have 
negative values. For isotropic elements = 2 XT^* XC 

anisotropic elements Fj ^2 input by the user. For composite materials it is 

recommended to use Fj^ = 0 [6] . 

(c) Shear panel element E44: 

SR=^ (3.6-3) 

D_ _ 


0^2 is the computed shear stress in the element and 

allowable. For stress constraints the constant function is defined as 


g = 1 - SR > 0 


(3.6-4) 


Fully Stressed Design (FSD) is applicable for the following elements: 

(a) Axial bar elements E23 

(b) Two dimensional membrane elements E31, E41 

(c) Shear panel elements E44 

For the bar and shear panel elements the stress ratios are the resize 
factors , but for the membrane elements the resize factors are calculated 
as follows : 
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B = (-ij- t 

^ XT XC YT YC 


-(o^)‘ 


C. = . 

^ XT^ • XC^ 


(oif ai ^ 

~r — r <■?> 

YT • YC S^2 


* 21^2 <’1 • 4 


Resize factor for the i~th elements, is 


(3.6-5) 


(3.6-6) 


12 ^''2 

RFi = I (B^ + (Bf + 1+C^) ) 


(3.6-T) 


When TFSD is applied the resize factors are calculated as follows: 
(a) Axial "bar elements E23: 


RF = 


M 


^1 -°t 


(3.6-8) 


where a., is the mechanical stress in the element. O. is the 

thermal stress and is either the tensile or compressive 

allowable depending on the sign of 
(b) Two dimensional membrane elements E31 


^XT XC^*^lt ^YT YC^*^2t 


(3.6-9) 


wKm ^ ^ y^2M 


(3.6-10) 
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2 


CT = - ---g*.- + (^) + 

XT ' XC YT • YC ■' 


12 


+ 2*F *0 *0 

12 It 2t 


(3.6-11) 


CM = - - 7r ~~~ + + 


XT • XC YT • YC 


12 


+ 2F O *0 
•^12 IM 2M 


( 3 . 6 - 12 ) 


Per GT Pa • a Pa • a 
__ It IM 2t °2M . 12t 12M . 

DMT = — - — _ — : — + ::: + 


XT'XC 


YT«YC 


(Si,)‘ 


+ 2F (o *0 +0 ' O ) 

12^ IM 2t 2M It^ 


(3.6-13) 


A = (BM + DMT)^ " It • CM (BT + CT - l) 


(3.6-llt) 


RF = J-(_(BM + DMT) + (A)^^^]/2(BT + CT - l) (3.6-15) 

1,2 2 - 

and the resize factor is the minimum nonnegative root. If both RF^ 

& RF^ are negative, the program sets RF to 2 . 

If the denominator of Eq. 15 is zero. 


RF 


- CM 
BM + DMT 


(3.6-16) 


in Eqs . 3.6-9 to 3.6-13 subscript t and M are for thermal and mechanical stress 

respectively. 
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(c) Shear panel elements EUU: 

RF = - (3.6-17) 

^12 ■ ‘^lat 

where o..,, & are the mechanical and thermal shear stress in the 

Ic-J"! 12t 

element . 

When one of the resize factors in the bar or shear panel elements 
is negative, the stress ratio is used for resizing this element. 

Input Data Sets; 

1. STRS EXY iset icase; Generated by GSF in library QLIB and contains 
the stresses and stress resxiltants. 

2. GOST EXY iset icase: Generated by DGRA in library GSLIB & contains 

the derivatives of the stresses. This data set is read when 
derivatives of stress ratios are required 

3. Secname BTAB MASK ID: Generated by TAB in library BLIB and 

contains section properties. 

U. ELTS ISCT 0 0: Generated by TAB in library BLIB and contains section 

identification number of all the element types in the problem. 

5. ISCT EXY 0 ID: Generated by KSDM in library DLIB and contains 

the section niimbers that control the elements. 

6. LOAD ATflT 00: To be generated by the user (by AUS) in library ALST 
when there is more than one set of allowable stresses. NI = number 

of load cases, NJ = 1. Contains the allowable set number corresponding 
to each load case. When FSD or TFSD resizing is applied, SCOII reads 

also the following data sets: 


38 



7 . 


KSDV Secnane 0 ID: Generated by DESVAR in library DLIB and contains 

the design variable numbers which control every section. FSD and 
TFSD assume that every section is controlled by a single design 
variable. 

8. FSDV Secname 0 ID: Generated by DESVAR in library DLIB and contains 

the corresponding factors of the design variables , see DESVAR for 
details . 

9 . DESI VARI 0 0: Generated by DESVAR or AUS in library DLIB and 

contains the current values of the design variables . 

10. CDIM Secname 0 ID: Generated by KSDM in library DLIB and contains 

the current structural dimensions of the sections. 

11. GDIM Secname 0 ID; Generated by KSDM in library DLIB and contains , 
the minimum gage values of the structural dimensions . 

12. MASS DER 0 0: Generated by KSDM in library DLIB and contains mass 

derivatives of the structural mass with respect to the design variables 

13. MING BW 0 0: Generated by KSDM in library DLIB and contain the 

minimum gage structural weight . 

To apply TFSD rather than FSD the user has to generate the data set 
MCTR LCAS 0 0 data set in library QLIB by using AUS. (NI = number of 
load cases, NJ = l). Each thermal and mechanical case should have a 

load case number and this data set contains for each mechanical 
load case the corresponding thermal load case, 
corresponds to purely mechanical load case 
0 correspond to thermal load case. 
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Output Data Sets: 

1. EQST EXY iset icase: Generated in library QLIB and contains the 

magnitude of the eqxii valent stress ratio defined earlier at the 
centroid of each element in a given element type EXY. For the two 
dimensionel elements this data set contains the six coefficients of 
Tsai Wu theory (Eq. 3 . 6 - 2 ). 

2 . If the allowable stresses are given as input in the input card stream 
ALST Secname set icase data set is generated as output in library 
ALST. 

3. RESZ EXY iset lease: Generated in library RLIB when FSD or 

TFSD is applied and contains .resize factors of the elements. 

1*. GGl CON set idcl: Generated into library QLIB and contains the 

critical constraints (g ^GLIM ) for each load case. The last name 
of the data set is equal to the load ceise n\miber plus the reset 
control COID. This data set is generated only if LFSD = 0. 

5. DGl CON set idcl: Generated into library GSLIB and contains the 

derivatives of the critical constraints. 

When FSD or TFSD is applied data set DESl VARl 0 0 is updated 

according to the resize factors or stress ratios calculated by SCON. 

Reset Controls : 


Name 

Default Value 

Description 

ELIB 

1 

Library containing element information 



packets . 

QLIB 

1 

Source library for displacement and 



stress data sets. 


ho 


Eeset 

Naae 

11 


Controls (Continued): 
Defatlt Veilue 


L2 


MSI 

CON 

G3LI3 

5115 


LF3D 


1 

1 

1 

1 


Description 

The first load case nvunber. In default 
it assumes the lowest load case niunber 
available for displacement data set. 

The last load case nvnnber. In default 
it assimes the highest load case for 
which displacement data are generated 
earlier. 

Load set number. 

Library containing the allowable 
stress data. 

Constraint case number. 

So\irce library for derivatives of 
stress data set. 

Source library for TAB and ELD 
generated data. 

Library containing the data sets 
created by DESVAR and KSDM. 

= 0 calcvilate stresses constraints 
and their derivatives 
= 1 calculate only stress ratios 
= 2 calculate stress ratios and update 
design variables using PSD or TFSD* 


*TPSD is used when the data set MCTR LCAS 0 0 exists in library QLIB. 



Reset Controls (Continued): 


Name 

RLIB 

GLIM 


Defa\ilt Value 
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COID 


Description 

Library containing resize factors 
calculated by FSD or TFSD 
Critical constraint cutoff point 
only constraints satisfying g ^ GLIM 
are stored in data set GGl CON 
Constraint set identification. 

Last name (NU) of data set GGl CON 
set idcl is equal to COID plus load 
case number. No other constraint 
set may have same last name. 


Execution Sequence: 

The following is the sequence of data cards 
@XQT SCON 

Reset Cards (if any) 

Stress allowables input as follows: 


STRESS 


Section type^^\ Set number 

ns^, ms^, XT^^\ XC^^\ 
(2) (2) 

nSg 5 mSg 9 9 9 

S^2^ YC^^^ 

s(2)^ YT^^\ 


TIC me 

nSj^ ^ niSj^ 5 XT ) XC 9 ^12 * 
Section type^^\ Set number^^^ 


Section type^^\ Set nxanber^^^ 


k2 


If there is only one set of sdlowable stresses y the set number is 
not required. But for solving a problem with more than one set, a data 
set LOAD ALST 0 0 should be generated by AUS in library ALST, (see Input 
Data Sets section). 

The meaning of the line: 
ns, ms, XT, XC, XT, YC 

is that sections ns to ms have the following allowable stresses; 

XT Allowable tensile stress in X direction 

XC Allowable compressive stress (negative value) in X direction 
Allowable shear stress 

YT Allowable tensile stress in Y direction 

YC Allowable compressive stress (negative value) in Y direction 

YT & YC are used only for the two dimensional elements E31 E32 E33 

E41 E42 e43. 

If the allowable stresses correspond to a single section ms is to be 
omitted. 
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SECTION U 

FLUTTER ANALYSIS AND DERIVATIVES 


Five processors AERO, SUBKRN, GAF, DKM and FLT are used to calculate 
of aerodynamic data, to perform flutter analyses and for obtaining 
derivatives of flutter constraints. 

Processors AERO, SUBKRN and GAF define the aerodynamic parameters 
and produce the generalized aerodynamic forces. They should always be 
executed in the order shown above. FLT performs the flutter analysis 
(I0PT=0 or IOPT=l) and (if I0PT=2) also calculates the derivatives of 
the flutter constraints with respect to the design variables . The 
processor DKM, which calcTilates derivatives of the stiffness and mass 
matrices with respect to the design variables, is retjuired only if 
^ 0 j"ivatives of flutter constraints are to be calculated in FLT. 

typical seq,uences in which these processors are used are 
given below. 


a. Flutter analysis 
TAB 


SPAR processors needed to obtain vibration modes 


EIG 
AERO 
SUBKRN 
GAF 
FLT 

b. Flutter anaJ.ysis and one resizing iteration 


TAB 


SPAR processors needed to obtain vibration 


modes 


DESVAR 


define design variables 
KSDM define minimum gages (with RESET parameter MING = O) 

AERO 

SUBKRN 

GAF 

DKM 

FLT 

NEWTON calculate search direction and' new values of design 

# 

variables 

KSDM calciilate structural sizes corresponding to new 

design variables 

The aerodynamic processors AERO, SUBKRN and GAF and the flutter 
processor FLT are based to a large extent on the WIDOWAC program (Ref. 2, 

7 ). 

U.l AERO - AERODTOAMIC INPUT PROCESSOR 

Function : AERO reads and processes the planform and aerodynamic grid 

data associated with a wing structure. The wing is assumed to be in the x, 
y plane, the chordwise coordinate is x and the spanwise coordinate is y. 

The aerodynamic grid consists of integration points vised for 
calculating the generalized aerodynamic forces (GAF) and collocation 
points used for solving the pressure-downwash integral equation (see 
Ref. 1). These points lie along chordwise lines whose number (NYS) and 
location are specified by the user. All the chords are used for 
integration but only a subset is used for collocation. The user has 
to specify the number (NYC) and identity of the collocation chords and 
the number of collocation points per chord (NS). Usually, both NS 
and NYC are either four or six. When optimization is to be performed 



it is reconjmended that the effect of these ninnhers (NS and NYC) on the 
flutter results he checked and the lowest values which yield acceptable 
accuracy he chosen because the optimization time is quite sensitive 
to these numbers. 

The vibration modes which are used in the ceilculation of the 
generalized aerodynamic forces are available at the Joints of the 
finite element structural model. A two dimensional spline is used to 
interpolate between the finite element nodes and the collocation and 
integration points. The user specifies the list of structural joints 
that are to be used to construct the spline for each chord. One of the 
considerations in choosing chord locations for the aerodynamic grid is 
the proximity of enough structural joints so that interpolation errors 
are minimized. If all the structural nodes lie on one chord the 
processor artifically doubles their number by constructing an image 
set at an adjacent spanwise station. For this reason NXP (see RESET 
'controls) must be doubled in this case. 

Reset Controls 

Name Default Value Meaning 

ALiB 1 Destination library for 

aerodynamic data. 

BLIB 1 Library containing the structure 

definition sets generated by TAB 

IHj ■ 0 Number of leading edge cranks* 

(Maximum - eight) 

NT 0 Nianber of trailing edge cranks* 

(Maximum - eight) 

■*^See discussion of namelist WNGEOM 


h6 


Name 

Defavilt Value 

Meaning 

NYS 

7 

Number of input spanwise stations+ 

NYC 

1* 

Number of spanwise stations used for 
collocation (less than or equal to 
NYS)+ 

NS 

1* 

Number of collocation points per 
chord+ 

NXP 

10 

Maxim\mi number of structural joints 
per spanwise station (note that if all 
nodes have same y coordinate NXP must 
be doubled) 


The reset controls are followed hy the following sets of data: 

1. PLAN - Planform definition 

A card consisting of the word PLAl^ is followed by namelist WITGEOM 

Namelist WNGEOM 


This namelist defines 

XLR 

XLT 

XTR 

XTT 

YTIP* 

XL ( ) 


the wing planform; 

X coordinate , leading-edge root 
X coordinate, leading-edge tip 
X coordinate, trailing-edge root 
X coordinate , trailing-edge tip 
y coordinate , tip 

a vector of length NL containing the x coordinates of 


crank points on leading edge 

YL ( ) a vector of length NL containing the y coordinates of 

crank points on leading edge 
+ See discussion of aerodynamic grid 

* It is assumed that at the root y = 0 



XT ( ) a vector of length NT containing the x coordinates of 

crank points on trailing edge 

YT ( ) a v^ector of length NT containing the y coordinates of 

crank points on trailing edge 

For an uncranked wing, the following parameters may be xised in namelist 
WHGEOM instead of the preceding list: 

2 

AR aspect ratio. Span /Area 

area surface area, (Span) (Root chord + Tip Chord)/2 

ANGLEl leading-edge sweep angle, degrees 

TAPER taper ratio. Tip chord/Root chord 

2. AERO - Aerodynamic grid definition 

The purpose of this input data is to define the locations of the 
aerodynamic chords to be tised in the aerodynamic calculation and the set 
of structural Joints that lie in the vicinity of these chords. Five to 
ten such chords are usually adeq^uate for the aerodynamic calculations. 

Their location should be chosen to be uniformly distributed over the 
wing and to minimize the distance between the chordline and the 
structural joints that are used with it. 

A record consisting of the word AERO is followed by the following 
records : 

(i) Chord definition coordinates; NYS non-dimensional y coordinates 
(fractions of semispan), must be monotonically increasing 

(ii) Definition of collocation chords - NYC integers. It is 

recommended that even nvimbers of collocation chords be used. 

(iii) Node definition records, one per chord - the nimibers of structural 
Joints lying in the vicinity of the aerodynamic chords. 
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Input /Output Data Sets 


The processor AERO needs basic data sets generated by TAB to define 
the joint coordinates. It generates the following data sets containing 
the definition of the wing geometry and the aerodynamic grid. These are 
COL WING, BL76 WING, BLK 5 WING, SUBK CRNK, and AERO GRID 
INPUT EXAMPLE ; 

RESET ALIB= 3 , NXP=lU, NL= 1 , NT=1$ 

PLAN$ WING PLANPORM DATA FOLLOWING ON NAMELIST WNGEOM 
$WNGEOM XLR= 0 ., XLT= 1030 . ,XTR=ll 80 . , XTT=ll 80 . , 

YTIP=850. ,XL( 1 )= 515 . ,YL( 1 )= 1 * 25 . ,XT( 1 )= 1180 . ,YT(i)=H 25 .$END 

AERO$ AERODYNAMIC GRID DEFINITION 

0. .25 . 1+5 .61 .71+ .85 .93$ NONDIMENSIONAL CHORD 

STATIONS 

21+67$ SECOND, FOURTH, SIXTH, AND SEVENTH CHORD 
USED FOR COLLOCATION 

1231+567$ NODES 1 THRU 7 DEFINE THE FIRST CHORD 

10 11 12 13 II+ 15 

18 19 20 21 

22 23 2I+ 25 : 26 27 28 29$ 

30 31 32 33 : 3I+ 35 36 $ 

SUBKRN - SUBSONIC KERNEL FUNCTION AERODYNAMICS PROCESSOR 
Fimction 

SUBKRN calculates aerodynamic matrices based on a subsonic kernel 
function formvilation. 
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RESET Controls 


Name 

Default Valxae 

Meaning 

ALIB 

1 

Soxirce and destination library for 

aerodynamic data reduced-frequencies 
set number 

The reset controls are followed by a card containing the word KSET and 

namelist NAMI 
The namelist 

Namelist NAMI 

defines the reduced frequency set 

Name 

Default Value 

Meaning 

MACH 

0. 

Mach nimiber 

KMAX 

U. 

Maximal reduced frequency 

KMIN 

0. 

Minimal reduced frequency 

KMED 

1. 

Median reduced frequency 

VK ( ) 
NK 

20 

Values of reduced frequencies 
Number of reduced frequencies 


The reduced frequency is defined as k = (h w/V) where h is the root semi- 

o o 

chord of the wing, u is a frequency , and V is the air speed. There are 
two options for defining the reduced frequency set. One option is to read 
in the values as the vector VK. If no values axe input into VK SUBKRN 
calculates the aerodynamic matrices for a set of NK reduced frequencies, 
half of which fall between KMIN and KMED and half between KMED and KMAX. 
There is no need to calculate the aerodynamic matrices for a large number 
of reduced matrices. Rather, it is. more economical to have NK on the 
order of 5-20, eind then \ise interpolation (in the flutter analysis) to 
obtain values for intermediate values of the reduced frequencies. 
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Inptit /Output Data Sets 

Processor SUBKRN generates two data sets using only input data sets produced 
by processor AERO: 

1) KSET bbbb* kset 0 

is of length 2*NK and contains the vector VK of the reduced 
freq.uencies and a vector VDK of the same frequencies mapped into 
the interval (-1,1). The variable kset is the frequencies set 
number 

2) AERO MATR kset 0 

contains mode independent aerodynamic matrices in NK blocks. 

Each block contains a real and imaginary matrix of order NS*NYC 
(see processor AERO), that is, the total number of collocation points. 

1|.3 GAP - GENERALIZED AERODYNAI-HC FORCE PROCESSOR 
Fvinction 

GAF calciilates the generalized aerodynamic forces and their derivatives 
with respect to the reduced frequencies for a given set of reduced 
frequencies. 

The generalized aerodynamic forces calculated in GAF are the array 

A. , = C ff Ap*^ (x,y) z^ds 
S 

where S is the wing area, Ap"^ is the pressvure differential associated 
with the j-th vibration mode, and z^ is the z component of the i-th mode 
(the wing middle— surface is assumed to be* close to the x— y place). The 
constant C 

* bbbb denotes a blank name 
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C = 8tt ql^/b 
^ o 

where q Is the dynamic pressure, 1 is the wing semi-span and the 

rcxjt semi -chord. 

Reset Controls 

Name Default Value Meaning 

Source library for data sets generated 
by AERO and SUBKRN 
Source library for vibration modes 
Number of vibration modes to be used from 
data set VIBR MODE 
Reduced frequencies set number 
First vibration mode to be used from 
VIBR MODE (that is, modes LI, Ll+1,. . . 
are used). 

BLIB 1 Source library for data sets generated 

in TAB 

Input /Output Data Sets 

Processor GAF uses most of the data sets generated in processors AERO and 
SUBKRN. Also, data sets defining the joint coordinates and active degrees 
of freedom from TAB and the vibration modes data set VIBR mode are 
needed. 

The only data set generated by GAF is GAF bbbb* kset 0. This data set 
contains two blocks of real and imaginary generalized aerodynamic forces, 
respectively. These have the dimension (NK, NEIG, NEIG) where NK is the 
number of reduced frequencies. This data set is written back into 
library ALIB. 

* bbbb denotes a blank name 


ALIB 1 

VLIB 1 

NEIG 2 

KSET 1 

LI 1 
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h.h DKM - GENERALIZED STIFFNESS AND MASS MATRIX DERIVATIVE PROCESSOR 
Fun c cion 

This processor is used to calculate derivatives of the generalized 
stiffness and mass matrices with respect to the design variables. 

Reset Controls 


Name 

Default Value 

Meaning 

ELIB 

1 

source library for element packets 

VLIB 

1 

source library for vibration modes 

DLIB 

1 

source library for design variables 

DKML 

1 

destination library for derivatives of 
generalized stiffhess and mass matrix 
blocks 

LI 

1 

First mode to be losed from data set 
VIBR MODE MASK MASK 

NEIG 

2 

Number of vibration modes data 

MXGR 

100 

Maximijm number of section for any section 
type 

MXGD 

5 

Maximum number of design variables that 
control the thickness of the same section 

NC 

2 

Nimiber of design variables 

GDIV 

1 

GDIV should be set to the value of the 

2 

acceleration of gravity (e.g. 386. U in/s 
if weight rather than mass data was 
supplied in TAB. Otherwise GDIV=1. 

DKM can 

handle only element types 

E23, E31, Eltl, and EUU. 
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Input /Oiit'ptit Data Sets 

The processor DKM reads most of the data sets generated by processor 
DESVAR. Additionally it requires data sets generated by ELD, E and EKS 
defining element properties and the vibration mode data set VIBR MODE. 

Two data sets are generated by DKM, 

1) DKM bbbb* 0 0 

containing the array DKM (NEIG, HEIG,NC+l). The first NC matrices stored 
in DKM are the derivatives of the generalized stiffness matrix with 
respect to the design variable. The last matrix is the generalized 
stiffness matrix obtained when all design variables are set to zero. 

2) DGM bbbb* 00 

This data set contains derivatives of the generalized mass matrix with 
respect to the design variables and has the same structure as the DKM 
data set . 

U.5 FLT - FLUTTER ANALYSIS PROCESSOR 
Fimction 

This processor may be used to perform flutter analyses and 
to calculate derivatives of flutter constraint functions. The vibration 
modes of the structure are used as generalized coordinates for the 
flutter analysis. When the structxire is resized new modes are calculated 
as a linear combination of the modes of the initial structure. For this 
process to be effective the number of nodes calculated for the initial 
structure should be considerably larger than the number of modes needed for 
* bbbb denotes a blank name 
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the flutter analysis. Typical numbers are 15 to 20 modes for the initial 
structure (KEIG) and 5 to 8 modes for the flutter analysis (NM). 

The flutter constraint used in the optimization process is the 
continuous flutter constraint of Ref. 3. This constraint formulation 
includes not only the flutter point but also "hump” mode critical points 
as constraints. The critical points on the V— g diagram are referred to 
in the output as "danger points". 

Reset Controls 

Meaning 

Source library for aerodynamic data 
Source library for design variables 
Source library for vibration modes 
Source library for derivatives of mass and 
stiffness matrices 
Riomber of design variables 
Three times the maximum number of danger 
points on the V-g diagram per flutter 
condition* 

Number of vibration modes for initial 
structure (order of GAF and DKM matrices) 
Number of modes needed for flutter analysis 
First vibration mode used in flutter 
analysis 

qLIB 1 Destination library for flutter constraints 


Name 

ALIB 

DLIB 

VLIB 

DKML 

NC 

MXDP 


NEIG 

NM 

LI 


Defavilt Value 
1 
1 
1 
1 

0 

. 6 


2 

1 


* For each danger point on the V-g diagram the two closest points are also 
included as danger points. 
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Reset Controls (Continued) 
Name Defaiilt Value 

HALF 1 


NF 

1 

SHIFT 

0 . 


lOPT 


1 


LUKIT 1 


MUHIT 


1 


CASE 


1 


Meaning 

and their derivatives 
Set to 0.5 if only one half of one ving 
structure is modeled (aerodynamic forces 
are multiplied by HALF) 

Niomber of flutter conditions 
Eigenvalue shift, used when rigid body 
modes are present; usually the sane 
value as used in EIG 

Analysis /optimization flag. I0PT=1 for 
V-g analysis starting from initial altitude 
and moving toward matched altitude. 

IOPT=0 for V-g analysis at a single altitude 
I0PT=2 same as I0PT=0 plus calculation of 
derivatives of flutter constraints. 

Ratio of user length unit to one meter. 

For example, if the length imit is one 
foot LUNIT=0.30it8 

Ratio of user mass unit to 1 kilogram. 

For example, if the mass unit is the 
slug, MUNIT=1U.6 

Last name of flutter constraint data set - 


GGl cor 1 CASE 

The reset controls are followed by a card containing the word FLUTTER 
and then namelist FLUTER which is repeated IlF times. 
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Namelist FLUTER 


Name 

H 

REF 

IGAIN 


Defavilt Value 


Meaning 

Initial (for analysis) constraint’ (for 
optimization) altitude 


1. rad/sec. Reference frequency 


10 


KMED 


1 . 


SG 

KSET 

MACH 


0.02 

1 . 

0 . 


Ratio of number of entries in V-g table 
to NK (number of reduced frequencies 
for which GAF were con 5 >uted) 

Median reduced frequency (needed only if 
automatic reduced frequency generation 
was called for in processor SUBKRN) 
Structural damping coefficient. 

Reduced frequency set to be used in 
flutter analyses 
Mach number 


Namelist FLUTER is repeated NF tines and then followed by namelist 
OPTIMUM for design variables values and increments. 

Namelist OPTIMUM 

Name Defaxilt Value Meaning 

0. Vector of design variables* 


C(I) 

DC(I) 


10 


-It 


Vector of design variable increments used 
for finite difference calculations of 
derivatives. On UNIVAC DC (l) shotild be 
about 1000 times smaller than the magnitude 
of the design variables. On CDC-10^ 
times smaller. 

* Required only if user wishes to override values obtained from data set 
DESI VAEI 00 
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InPTib /Output Data Sets 

Processor FLT vises most of the data sets generated by processors AERO, 
SUBKRN, GAP and DKM. Additionally the values of the design variables, 
DESI VARI and the vibration frequencies VIBR EVAL MASK MASK are used. 

The data sets GGl CON 1 case and DGl CON 1 case containing, respectively, 
the values of the flutter constraints (as defined in Ref. 7) and their 
derivatives with respect to the design variables are generated by FLT. 
Two data sets containing the definition of the flutter conditions, FLTR 
PACK 0 0 and FLTR PACK -1 0, are also generated and written back into 
library ALIB, 
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SECTION 5 
OPTIMIZATION 


The PARS processors give the user the flexibility of having either 
automated optimization or "a man in the loop” process. The automated 
optimization is furnished through the use of loop processors, SOL and FOL 
and FSDL that call other processors in a predetermined sequence. The 
static optimization loop SOL and the flutter optimization loop FOL use 
mathematical programming for the optimization and call processor NEWTON 
for the resizing process. The third loop, FSDL, provides the user with 
the option of fully stressed design resizing when only stress constraints 
are considered. The viser has the option of using these optimization 
loops or creating his own loops by calling the processors in sequence. 

The sequence does not have to be the same as the one provided by SOL, FSDL 
or FOL. The following are two examples of such nonstandard sequences. 

a. Simultaneous optimization for some outside generated constraints 
say for example buckling, (which passes the information about the 
constraints and its derivatives throu^ data sets) and some internally 
generated constraints like stress, displacement or thermal constraints. 

It is assumed that BUCK is the name of the group of some outside 
processors which generate a set of buckling constraints and its 
derivatives during each new analysis process. We can create the 
following loop which performs the simultaneous optimization. 

@XQT TAB 

, Static analysis 

§XQT SSOL 59 



§XQT GSF 
XQT PSF 

§XQT DESVAR $ Define design variables 
§XQT KSDM $ Define minimiua gage 
* @XQT BUCK $ BUCKLING Analysis and derivatives 
@XQT DGRAD 

@XQT SGRAD Displacement and stress constraints 
§XQT DCON and their derivatives 
@XQT NEWTON $ Resize 
@XQT KSDM $ Update structural sizes 
@XQT E 
• 

@XQT SSOL 
@XQT GSF 

Repeat sequence from (*) 

b. Several one-dimensional searches based on one analysis. 

The user may elect (following t' ' '^SS code [8]) to resize the 
structure several times based on ' s of the constraints 

obtained from processors DCON SCON and BUCK rather than repeat the 
analysis and derivative calculations after each one dimensional search as 
is done in processor SOL. This only requires calling NEWTON several 
times in succession. 
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5.1 NEWTON : Newton’s Method Optimizer 

Function: NEWTON finds a search direction and the minimvim in that direction 

in a sequence of xmconstrained minimizations (SIJMT) , vising Newton's -method. 

An extended interior penalty function, is employed to account for constraints. 
Even though NEWTON performs only a single one dimensional search it 
performs convergence tests which are used when this processor is used in a 
loop. 

• V 

Optimization procedure J ' 

The object is to minimizeHhe mass M of the structvire subject to 
the constraints 

g^(v) ^0 (i, . . .maxcon) (5.1-1) 


where , 

V represents the vector of n design variables, g^ is the ith 
constraint and maxcon is the number of constraints. The constrained 
optimization problem is transformed into the following series of 
unconstrained problems: 

Minimize 

i=maxcon 

P(r, v) = M (v) + r E f. (v) - (5.1-2) 

i=l 


for r^ , Tg , 


-»■ 0 


The sm which is added to the mass is called the penalty function. 

The definition of f^ in (5.1-2) for a quadratic extended penalty function is 
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(5.1-3) 


if Si > 6^ 


(v) =( 


- [(g./g - 3(g,/g^) + 3] if g, < g 




1 O 


1 — O 


where g is a transition point* 
o 

The definition of f^ in (5*1-2) for a cubic extended interior penalty 

function [9] formulation is 


ri 


fl (V) . 


if Si>8„ 
2 


{5.1-h) 


- [A (- 1 + — ) + (1-^) +)1-^) + 1] 

vSq So ®o So 


if gi Ig,. 


The transition parameter changes with r as 


So = S„i(r/r^) 


where and are the initial values of g^ and r, respectively, and 

p is a constant (usually 0.5). For each value of r the unconstrained 
minimization is performed using approximate second derivatives reqviired 
for generating a suitable search direction. A one dimensional search is 
then conducted for locating the minimum in that direction. The minimization 
for- a single value of r is terminated when any one of the following 
criteria is satisfied. 
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(1) The Euclidean norm of the gradient of P is smaller than 0.05 
times the gradient of M. (gradient convergence). 

(2) NEWTON estimates that at the beginning of the last one- ■ 
dimensional search P(r,v) is within 2 % from its minimum value (Quadratic 
convergence ) . 

(3) LCNT (a value controlled by a reset parameter) number of 
one— dimensional searches have been executed. 

When one of the convergence criteria is met r is reduced by a factor 
RDC (see RESET controls) unless both of the following convergence criteria 

are met 

(1) The penalty fianction is less than one percent of the weight 

(2) The reduction in the weight during the last one dimensional search 
is smaller than one percent . 

Constraint Approximation 

NE\^T0N performs a complete one dimensional search without calling any 

analysis routines by approximating a constraint function g(v) based on 

their values at the initial point of the search v^ (with components 

Y i=l,...n) and their derivatives at that point. The approximation is 
oi 

(5.1-5) 

(5.1-6) 


g(v) s E(v^) * E If- 

J 


where 


if ag/av. < 0 
J 


B. = 
J 


V . + V . 

m.i 0,1 

V , + V 

mj j 


if ag/av. > 0 

j 
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and V . is a vector of constants input by the user (as GAGE in processor 
mj 

DESVAR). Hie value B = 1. corresponds to a linear approximation and the 

J 

alternative expression in Eq. (5* 1-6) corresponds to an inverse approximation 
The choice is based on which one is more conservative (see [10], that is, 
which choice results in a smaller value of g(v). Because in PARS a zero 
value for the design variable usually corresponds to a nonzero value of the 


structural parameter that it controls (see Eq. 2,1-1) the constants V 


mj 


are introduced. In the simple case where a single design variable 
controls a structural parameter q^ , Eq. (2.1-1) becomes 


q. = q . + V.F. . 
^oj 1 ij 


In such a case would correspond to an approximation 

linear in the reciprocal of q.. In the more general case where a 

J 

single design variable controls several structural parameters or vice 


versa the choice of a good value of is more difficult. The guideline 

to that choice is that when Vj = 'the approximation becomes unbounded. 

For example if 'V^ controls the thickness of one lamina in a composite 
material a reasonable choice for 'Vj^j is the laminate thickness when 
■Vj = 0. If a very large value is selected for the approximation 

becomes linear. 


Relation to Other Processors 

NEWTON can be used as an integral part of PARS or as a stand alone 
processor used in conjunction with another application progrEim that 
calculates the constraints and their derivatives . If NEWTON is used in 
such a mode the user, through the use of processor AUS or through the 
application program, needs to supply data sets which are otherwise 
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generated by other PARS processors. Ihe folloving data sets ara required 
by NEWTON. 

1. DESI VARI 0 0: Contains the values of the design variables at the 

beginning of the search. On exit it contains the updated value of the 

design variables. (NI = 1, NJ = NDV) 

2. DVMX VARI 0 0: Contains the maximum values of the design variables 

and is required only if such bounds are used. (NI = 1, NJ = NDV) 

3. GAGE VARI 00: Contains the values of V . (NI = 1, NJ = NDV) 

mi 

U. MASS DER 0 0: Contains the derivatives of the mass with respect to 

the design variables. 

5. GGl CON set case: Contains the values of the constraint functions 

g^ (NI = number of constraints, NJ = l). 

6. DGl CON set case contains the derivatives of the constraints with 
respect to the design variables (NI = NDV, NJ = 1, each constraint 
in a block) . 

NEWTON creates a data set D12A BLOC 0 0 which contains the parameters 
that control the optimization, most of which are reset parameters. If 
this data set exist it overrides the reset values. 

The contents of this data set are: 

R ,RDC , GLIM , GMIN ,P , GIC , DELC ,LCNT ,K0UNTR ,NFINAL where 

GMIN - the value of the most critical constraint function excluding design 
variable bounds 

KOUNTR - the number of the present one dimensional search in the sequence 
for the current r- value 
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HFINAL - a teraiinatiou flag. Set to one ty NEWl'ON vhen termination criteria 


are met. 

The other parameters in the data set axe explained as reset parameters. 


Name 


LCNT 


R 


RDC 

GLIM 

P 

GIC 

DC 


Default Values 


T* 


50 . 


20 .* 

0.5 

0 . 5 * 
0 . 32 * 

1. E-04 


Description 

LCNT is the maximum number of one 
dimensional searches 

Factor of penalty function [r]., recommended 
value, r=0.1 M where M is the expected 
minimtmi mass. 

Reduction factor for r [l/RDC]. 

Limiting value for critical constraints. 

Exponent of [ = cr J . 

The initial value for g . 

o 

The increment in finite difference 

api 


HHLI 


Destination library for Hession matrix. 


GGLI 


1 


SET+ 

L1+ 

L2+ 

DVLIB 


1 

1 

1 

1 


OEPF 


2 


Library containing GGl CON mask mask 
and DGl CON mask mask data sets. 
Constraint set. 

First case constraint. 

Last constraint case. 

Library containing data sets generated 
by processor DESVAR. 

Order of extended interior penalty 
function to be used. 


* Default value is recommended. 

+ Newtoti will process all constraints stored in data sets GGl CON SET CASE 
where CASE runs from LI through L2. 
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5 SOL - Static Optimization Loop 

Function; The SOL processor provides a structural optimization loop 
for static-strength (stress and displacements) and minimum-gage constraints 
under several alternative load conditions. Various processors (the list is 
given below) are called in a sequence to obtain the minimum weight structural 
design. The user may elect to rim the individual processors one by one 
instead of using SOL. This, however, will require repeating the reset • 
cards for each call. 

It is assumed that minimum gage and static stress constraints are 
always present. Constraints imposed on joint displacements may or may not 
be specified. Design variables can be a physical dimension of a 
structural element (thickness for 3 and h node elements and area for 2 node 
elements) or the nodal temperatures specified at certain joint locations. 

The following is the sequence in which various processors are called. 
Processor identification number and their functions are outlined briefly. 
Processor 

Processor Identification 

Name Number (IDP) Description 

KSDM 28 Updates structural sizes. 

E 35 Generates element packets 

EKS 36 Calculate intrinsic stresses and 

stiffness matrices 

K 13 Assembles system stiffness matrix. 

INV 16 Factors assembled system stiffness matrix. 

SSOL* 18 Computes displacements and reactions due 

to point loading applied at joints. 

* These processors correspond to standard SPAR processors , however, they are 
modified here to some extent in order to include them in the optimization loop. 
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Pi'ocessor 

Name 

Processor 
Identification 
Nimiber (IDP) 

Description 


EQNF* 

17 

Computes equivalent joint loading and 


GSF» 

20 

initial strain data. 

Generates element stresses and internal 


PSF« 

21 

load information. 

Prints element stresses and internal load 


DGRAD 

29 

information. 

Calculates derivatives of displacements 


TGRAD 

27 

with respect to design variables control- 
ling a thickness or area. 

Calculates derivatives of displacement with 


DCON 

32 

respect to design variables controlling 
nodal temperatures . 

Forms the displacement ratios for the 


SGRAD 

31 

nodes where nodal displacements are 
constrained. 

Evaluate the derivatives of element 


SCON 

30 

stresses. 

Forms the stress ratio at the centroid 


NEV7T0N 

34 

of finite elements 

Implements the sequence of vinconstrained 

■ M 


minimization technique using Newton's 
Method. 


* These processors correspond to standard SPAR processors , however, 
they are modified here to some extent in order to include them in 
the .optimization loop. 
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One or more of the above processors can be skipped to restart the 

static optimization loop in case a run terminates in the middle due to 

some \inforeseen circumstances like maximum time or einy input data error. 

For detail refer to "Restart Execution Mode." later in this section. 

Reset Controls: 

Name Defavilt Value Description 

ELIB 1 Library containing the element packets. 

QLIB 1 Library containing static displacement 

data sets. 

BLIB 1 Library containing output data sets 

generated by TAB and ELD processors. 

MIUG 0 (See processor KSDM) if a non zero 

value of KENG is specified in this 
processor the reset parameter HID 
should be set to a non zero value in 
processor KSDM. 

IPRT 0 If other than zero, processor PSF is 

called to print out element stresses and 
derivatives of element stresses. 

ITDV 2 ITDV=2 - Problem does not contain temperatures 

ITDV=1 - Temperatures are used but no 

temperature design variables 

ITDV=0 - Temperature design variable are used 
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Default Value 


Description 


AMAX 2 The maximum number of times the optimization 

loop is executed and the design variables 
are tqidated before the run is terminated. 

IDP 28 This is used to restart the loop if SOL 

terminates abnormally maximum time or 
otherwise. 

Execution Control Parameters: 

The following is the sequence of data cards required when executing 
SOL for the first time (for restart refer "Restart Execution Mode") 

@XQT SOL 

Reset Ceirds (if any) 

STOP$ IDF = 33, SOL 

V 

\ 

RESET IND=1 

a. (Additional Reset cards for KSDM) 

STOP$ IDP = 28 , KSDM 

b. (Reset cards for E) 

STOP $ IDP = 35 

c. (Reset cards for EKS) 

STOP $ IDP =36 

d. (Reset cards for K) 

STOP $ IDP = 13, K 

e. (Reset cards for INV) 

STOP$ IDP = 16 , INV 

f . (Reset cards for SSOL 

STOP$ IDP = 13, SSOT. 


TO 





(Reset cards for GSP) 

STOP$ IDP = 20, GSF 
(Reset cards for PSF) 

STOP$ IDP = 21, PSF 
(Reset cards for DGRAD) 

ST0P$ IDF = 29, DGRAD 
(Reset cards for TGRA) 

ST0P$ IDP = 27, TGRAD 
(Reset cards for DCON) 



These cards are optional and should 
be xxsed only if reset parameter 
»IPRT* 4 0. in ’SOL* 


DISCON: i^, adv^: ig, jg, adVgi i^, iy adv^ 

ST0P$ IDP = 32, DCON 


(Reset cards for SGRA) 

ST0P$ IDP = 31, SGRAD 
(Reset cards for DPSF) 

ST0P$ IDP = 23, DPSF 
(Reset cards for SCON) 

STRESS: EXY: ng^, mg^, s|^\ s^^\ s|g^ $ 

^^§ 2 ’ ^ 2 ’ ^ 12 ’* ^^ 3 ’ ^1 ’ ^2 ’ ^12 
$ 

ST0P$ IDP = 30, SCON 
(Reset cards for NEWTON) 

ST0P$ IDP = 3U, NEWTON 
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STOP$ card serves as an input data card separator for a processor. 
When it is encountered during a read mode, it stops scanning any further 
data cards. 

Restart Execution Mode: 

In restart , the reset and STOP$ cards for main processor SOL are 
always required, and should he supplied first. The other cards should 
follow the normal execution mode order. A few of the typical situations 
are illustrated through examples. 

Example 1. 

Assume that while executing DGRAD (item i of SOL) the run was 
terminated due to maximum page or time limits or may be due to improper 
RESET cards. If now one intends to restart the program, the following 
sequence shoixld be used: 

@XQT SOL 

RESET IDP = 29 

(Other Reset cards if any) 

STOP$ IDP = 33, SOL 

i. (Reset cards for DGRAD) 

STOP$ IDP = 29, DGRAD 

j. (Reset cards for TGRAD) 

ST0P$ IDP = 27, TGRAD 

k. (Reset cards for DCON) 

DISCON: i^; adv^: ig, ^2’ ^^''^2" ^3’ '^3* ^*^^3 

STOP$ IDP = 32, DCON 

l. (Reset cards for SGRAD) 


STOP$ IDP = 31, SGRAD 



m. 


n. 


} 


(Reset Ccirds for FSF) 

STOP$ IDF =23, PSF 
(Reset cards for SCON) 

STRESS: EXf: ng^, mgj_, 

( 2 ) ( 2 ) ( 2 ) 

’ ®2 ’ ^12 ’ “ 2 ^ 

STOP$ IDP = 30, SCON 
(Reset cards for NEWTON) 

ST0P$ IDP = 3U, NEWTON 


These cards are optional and should 
be i:ised only if reset parameter 'IPRT* 


i 0. in ’SOL' 

®2 ’ 12 ^ 

g(3) „(3) 


s^3)..„$ 
12 ^ 


Reset IND = 1 

a. ( additional Reset cards for KSDM) 

ST0P$ IDP = 28, KSDM 

b. (Reset cards for E) 

ST0P$ IDP = 35 

c. (Reset cards for EKS) 

ST0P$ IDP = 36 

d. (Reset cards for K) 

ST0P$ IDP = 13, K 

e. (Reset cards for INV) 

ST0P$ IDP = l6, INV 

f. (Reset cards for SSOL) 

ST0P$ IDP = l8, SSOL 


g* 


h. 


(Reset cards for GSF) 

STOP$ IDP = 20, GSP 
(Reset cards for PSF) 

ST0P$ IDP = 21, PSF 


These cards are optional and 
should be used only if reset 
parameter 'IPRT 0. in SOL. 


73 



Example 2. 

Assume that each one of the processors a to h has already heen 
executed separately outside SOL. To optimize the structure, the SOL 
loop can be entered at i using the same set of cards as shown for 
Exainple 1. 

5.3 FOL - FLUTTER OPTIMIZATION LOOP. 

FUNCTION ; 

The FOL processor is a flutter optimization loop for the minimvim weight 
design of a structure subject to flutter and minimvim gage constreiints . 

It calls the various processors required for the optimization and thus 
saves the tiser the tedium of calling them himself and repeating RESET and 
other input cards. Design variables are element thicknesses or areas 
(see DESVAR) and concentrated mass at structural joints. 

The following is the sequence in which various processors are called. 


Processor 

Name 

Processor 
ID Number 


Function 

KSDM 

28 

Updates intrinsic stiffness and stress 



matrices and global (diagonal) mass matrix 

E 

35 

Generates 

element packets 

EKS 

36 

Calculates 

intrinsic stiffness and stress 



matrices 


K 

13 

Assembles 

global stiffness matrix. 

AUS 

23 

Arithmetic 

Utility System* 

INV 

l6 

Factors assembled stiffhess matrix. 

EIG 

22 

Calculates 

vibration modes. 


* AUS is most commonly used for forming the desired mass matrix or 
introducing a shift into the stiffhess matrix. 



Processor 

Processor 


Name 

ID Number 

Function 

GAP 

38 

Calculates generalized aerodynamic forces . 

DKM 

uo 

Calculates derivatives of generalized 
stiffness and mass matrices. 

FLT 

39 

Performs flutter analysis and calculates 
derivatives of flutter constraints. 

NEWTON 

3U 

Updates design variables xjsing Newton's 



method with approximate second derivatives. 

RESET CONTROLS: 


Name 

Default Value 

Description 

BLIB 

1 

Library containing TAB & ELD generated data 
sets. 

ELIB 

1 

Library containing element packets. 

MLIB 

1 

Library containing mass matrix. 

MING 

0 

See processor (KSDM), if a nonzero value of 
MING is specified in this processor the 
reset parameter IND should be nonzero in KSIM 

IDP 

28 

ID of first processor called by FOL. 

AMAX 

2 

The number of times the processor loop is 
repeated. 

FRAC 

Oil 

Maximum fractional change in the square 


of the first vibration frequency that is 
acceptable without mode updating. 



VARIATION IM LOOP STRUCTURE; 


The "two reset parameters IDP and FRAC give the user the possibility 
of changing the basic loop structure. IDP is the identification number 
of the first processor to be executed and thus permits the user to start 
the loop at any point. This freedom is especially important for restart. 

FRAC permits the user to loop for most of the design process only over 
the last two processors FLT and NEWTON. The idea is that if the changes 
in the structure are not significant there is no need to update the vibration 
modes every time the design variables are changed. The change in the 
dynamic properties of the structure is measured by the change in the first 
vibration frequency. The user specifies FRAC which is the maximTim 
fractional change of the square of the first vibration frequency that may 
be tolerated without updating modes . 

EXECUTION CONTROL PARAMETERS: 

The following is the sequence of data cards required when executing 
FOL with IDP=28: 

@XQT FOL 
Reset Cards 
Stop$ FOL 

Reset Cards for KSDM 
Stop$ KSDM 
Reset Cards for E 
Stop$ E 

Reset cards for EKS 
Stop$ EKS 
Reset Cards for K 
Stop$ K 
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Reset Cards for AUS 


Input Cards for AUS 
Stop$ AUS 

Reset Cards for INV 
Stop$ INV 

Reset Cards for EIG 
Stop$ EIG 

Reset Cards for GAF 
Stop$ GAF 

Reset Cards for DKM 
Stop$ DKM 

Reset Cards for FLT 

Data for FLT (Namelist FLUTER, etc.) 

Stop$ FLT 

Reset Cards for NEWTON 
Stop$ NEWTON 

Repeat Input Cards for AUS 

Examples ; 

Several typical processor sequences are presented 
a.) Optimization of a wing structure which has heen analyzed in the past, 
@XQT TAB 

I 

Define structural grid 


@XQT ELD 
@XQT TOPO 
@XQT DESVAR $ 
@XQT AERO 
@XQT SUBKRN 
@XQT FOL 


Define design variahles 
Basic aerodynamic calc\ilation 
Flxitter optimization 
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b. Flutter Analysis and then optimization. 


gXQT TAB 
gXQT ELD 
gXQT TOPO 
gXQT E 
gXQT EKS 
gXQT K 
gXQT INV 
gXQT AUS 
gXQT AERO 
gXQT SUBKRII 
gXQT GAP 
gXQT FLT $ 
gXQT DESVAR $ 
gXQT SKDM $ 
gXQT FOL$ 


Define structviral grid 


Calculate vibration modes 


Calculate GAF's 

Flutter analysis 
Define design variables 

Define ninimim gage structure (set MING=0) 
Flutter optimization 


In this case the RESET cards following gXQT FOL are RESET IDP-UO, . . .$ 


loop starts at DK31 
stop$ FOL 

Reset cards for DK-I 
stop$ DKM 

Reset Cards for FLT 
Data for FLT 
stop$ FLT 

Reset cards for NEWTON • 

stop$ NEWTON 

Reset Cards for KSDM 

Btop$ KSDM 

• 

eards for GAF 7R 



c. Restart 


Suppose a previovis run stopped on maximum time after executing IRV, 
and we wish to restart. It is assumed that data sets have been preserved. 
@XQT FOL 

Reset cards for FOL (lDP=22) 
stop$ FOL 

Reset cards for EIG 

stop$ EIG 
• 

Reset cards for INV 
stop$ INV 

5.1* FSDL - Full Stress Design Loop 

The FSDL processor provides a Fully Stressed Design loop (FSD or TFSD) 
for static strength (applied forces and temperat'ures ) and minimum gage 
constraints under several alternative load cases . 

Various processors are called in a sequence to obtain the Fiilly 
Stress structural design. The \iser may elect to rxm the individual 
processors one by one instead of ■using FSDL. This, however, requires 
repeating the reset cards for each iteration. 

TFSD (Temperature Full Stress Design) algorithm, which often 
converges faster than FSD when thermal stresses are significant, is also 
available. If the data set MCTR LCAS 0 0 is foimd in QLIB library (see 
SCON processor for details) TFSD is used rather than FSD. 

The following is the sequence in which various processors are called: 
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Processor 

NflinR 

ISDM 

E 

EKS 

K 

INV 

EQNF 

SSOL 

GSF 

SCON 

Reset controls: 

Name 

MING 

AMAX 


IDP 


Processor 
Identification 
Number 

28 

35 

36 

13 

16 

17 

18 

20 

31 


Description 

Updates structural sizes.. 

Generates element rockets. 

Calculates intrinsic stresses and 
stiffness matrices. 

Assembles system stiffness matrix. 
Factors assembled system stiffness 
matrix. 

Computes equivalent joint loading 
and initial strains. 

Computes displacements and reactions 
due to applied joint loadings. 
Generates element stresses. 

Generates nondimensional stress 
ratios and resize factors. 


Default value Description 

0 See processor KSDM. 

2 Maximum of iterations FSD (or TFSD) 

loop is executed, design variables 
and ciurent section dimensions 
(CDIM data sets) are updated 
before termination. 

28 Used to restart the loop if FSDL 

terminates abnormally. 
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ITDV 


2 


= 1 EQKF is executed in the loop, 

= 2 EQNF is not executed in the loop, 
0 =0 Old data sets are being 

disabled and new ones (created 
by executing the different processors 
in the current iteration ) are 
generated in cumulative fashion. 

= 1 New data sets are generated in 
the same disk location as the 
old ones. 

The following is the sequence of data cards required when executing 
FSDL (on GDC computer) when terminal load cases are applied. 

gXQT FSDL 

Reset cards for FSDC 
ST0P$ FSDL IDP=33 
Reset IND=1 

Additional, reset cards for KSDM 

ST0P$ KSDM IDP=28 

Reset cards for E 

ST0P$ E IDP=35 

Reset cards for EKS 

ST0P$ EKS IDP=36 

Reset cards for K 

ST0P$ K IDP=13 

Reset cards for INV 

ST0P$ INV IDP=16 

Reset cards for SSOL 


INPL 

r 


1 
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STOP$ SSOL IDP=18 
Reset cards for GSF 


STOP$ GSF IDP=20 
Reset cards for SCON 

data for allowable stress (See in SCON processor) 
STOP$ SCON IDP=31 

When nodal temperatures are applied execute FSDL with reseting ITDV to 1, 
and add the reset control cards for EQNF and STOP$ card after the STOP$ 
card of INV. 
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